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H '. 1 Introduction 

"^ ■ Mass transfer in neutron star low-mass X-ray binaries (NS LMXBs) is either driven by 

I ■ loss of orbital angular momentum or nuclear evolution of the donor star, causing the 

^ . orbit to the shrink or expand respectively. The "bifurcation period" Pbif? the initial 

^ , binary orbital period which separates the formation of converging systems from the 

diverging systems (Tutukov et al. 1985) was found to be in the range ~ 0.4 — 0.7 day 

for LMXBs, and strongly dependent on magnetic braking (MB) (Pylyser & Savonije 

^ ■ 1988). For sufficiently small initial orbital periods, the evolution may lead to the 

0> , formation of of ultra-compact X-ray binaries (UCXBs) with orbital periods (P < 50 

QQ \ min), in which the donor is a white dwarf or a compact core of an evolved giant star 

(Nelson et al. 1986; Tutukov et al. 1987; Pylyser & Savonije 1988; Podsiadlowski et 

in ; al. 2002; van der Sluys et al. 2005). 

In this work we present the results on the evolution of NS LMXBs and the for- 
mation of UCXBs (Ma & Li 2009 for details). We consider the following processes 
related to mass and angular momentum loss mechanisms in LMXB evolution. (1) 
'k>( \ The standard MB law (Verbunt & Zwaan 1981; Rappaport et al. 1983) was shown to 

^ ; be contradicted with the observation of young stars in open clusters, and a modified 

version was proposed (Sills et al. 2000; Andronov et al. 2003). (2) There is strong 
evidence that during LMXB evolution the mass transfer is highly non-conservative. 
Possible ways of mass loss include "evaporation" of the donor (Ruderman et al. 1989) 
or "radio-ejection" of the transferred material (Burderi et al. 2001, 2002; D'Antona 
et al. 2006) due to the pulsar radiation/wind impinging on. In the latter case, the 
matter is lost from the system at the inner Lagrangian (Li) point, carrying away an- 
gular momentum and altering the period evolution. Additionally, a small fraction of 
the mass lost from the donor may form a circum-binary (CB) disk around the binary 
rather accretes onto the NS (van den Heuvel 1994). 



2 Evolution code and binary mode 

2.1 The stellar evolution code 

We use an updated version of the stellar evolution code originally developed by (Eggle- 
ton 1971, 1972) to calculate the evolutions of binaries consisting of an NS (of mass 
Ml) and an MS secondary (of mass M2). For the secondary star we assume a solar 
chemical composition {X = 0.70, Y = 0.28, and Z = 0.02). We assume that the 
spin of the secondary star and the binary orbital revolution are always synchronized. 
Assuming rigid body rotation of the secondary star and neglecting the spin angular 
momentum of the neutron star, the total angular momentum J of the binary system 
can be expressed as 

J = I2UJ + Jorb 

= l2io + G^/^MiM2iMi + M2y^^^'oj''^^^ (1) 

where I2 is the moment of inertia of the secondary star, G the gravitational constant, 
and u the angular velocity of the binary. 

We consider three kinds of mechanisms of angular momentum loss. The first is 
the angular momentum loss due to gravitational radiation 



dJcR _ 32 G^/^ MfMliMi + M 



1 1/2 



dt 5 c5 a7/2 



(2) 



where a is the orbital separation and c is the speed of light. 

The second angular momentum loss mechanism is for non-conservative mass trans- 
fer. We assume that a small fraction 6{<^ 1) of the mass lost from the donor feeds 
into the CB disk rather leaves the binary, which yields a mass injection rate of the 
CB disk as Mcb = —SM2. Tidal torques are then exerted on the binary by the CB 
disk via gravitational interaction, thus extracting the angular momentum from the 
binary system. The angular momentum loss rate via the CB disk is estimated to be 
(Taam & Spruit 2001) 

dJ, f2na^\ . / t \V3 

lcB = -7K^ Mcb - , (3) 



dt '"° V ^ / vt.. 

where 7^ = rja = 1.7 (r; is the inner radius of the CB disk), t is the time since mass 
transfer begins. In the standard a- viscosity disk (Shakura & Sunyaev 1973), the 
viscous timescale tvi at the inner edge r^ of the CB disk is given by t^i = 2'j^P/ 3710/3"^, 
where a is the viscosity parameter (we set a = 0.01 in the following calculations), 
/3 = Hi/ri ~ 0.03 (Belle et al. 2004), and Hi is the scale height of the disk. We also 
assume that the NS accretion rate is limited to the Eddington accretion rate, and 
that when the mass transfer rate is less than MEdd, half of the mass is accreted by 



the NS, i.e., Mi = min(MEddj — M2/2). The excess mass is lost in the vicinity of the 
NS through isotropic winds, carrying away the specific angular momentum of the NS, 
i.e. 

dJ r lM2alu, IM2I <2MEdd 



ML 



dt ''"'^ 1 (M2 + MEdd)a?w, IM2I > 2MEdd 



(4) 



where oi = aM2/{Mi + M2) is the orbital radius of the NS, and u is the orbital 
angular velocity of the binary. 

The third angular momentum loss mechanism is MB. We use the saturated mag- 
netic braking law suggested in Sills et al. (2000), 

^iMB-|-''"'(*)'^(ft)"''; '^-'^" (5) 

*l"=-|-A-..V(||)"=(||)-'^ .>.„ ''' 

where K = 2.7 x 10^^ gcm^s (Andronov et al. 2003), Wcr is the critical angular 
velocity at which the angular momentum loss rate reaches a saturated state, and can 
be estimated as ( (Krishnamurthi et al. 1997), 

Ucrit) =a;cr,0-^^, (6) 

n 

where Wcr,© = 2.9 x 10^^ Hz, rtg,© is the global turnover timescale for the convective 
envelope of the Sun at its current age, rt for the secondary at age t, solved by integrat- 
ing the inverse local convective velocity over the entire surface convective envelope 
(Kim & Demarque 1996). 

3 Results 

The formation and evolutionary paths of UCXBs depend on the adopted values of 
6. To illustrate the effects of S on the binary evolution, in Figs. 1 and 2, we plot 
the evolution of the donor mass and period as a function of age respectively, for a 
binary system with M2,i = I.IM©, Pj = 1.04 d and different values of 6. A larger 
value of S leads to shorter formation time, as seen from Fig. 2; if 6 is too small, the 
binary will not be able to reach the 50 min period within 13.7 Gyr due to inefficient 
angular momentum loss. When S < 0.0055, the orbital period first decreases with 
mass transfer until the donor star loses its outer envelope and shrinks rapidly at 
P ~ 0.1 — 0.2 d. This causes a cessation of mass transfer. In the subsequent evolution 
the orbital period may decrease down to the ultra-short regime under the effect of GR, 
until the secondary star fills its RL again, and the binary appears as a UCXB. When 
6 > 0.0055 the binary evolves directly into the ultra-short regime with decreasing 
orbital period. 
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Figure 1: Evolution of the donor mass for different values of the CB disk parameter 
5. 



We need to mention that the distribution of 6 depends on the value of the viscous 
parameter a. From Eqs. (3) one can see that the CB disk-induced angular momentum 
loss rate is proportional to a^^^5. So if keeping a^/^5 constant, the binary evolution 
will be exactly the same. 



4 Comparison with observations 

There are currently 10 UCXBs with known periods, 5 of which are persistent sources 
and 5 are transients. To compare observations with our CB disk-assisted binary 
model, we plot the Mi(= — M2/2) vs. Porb relations in Fig. ?? for binary systems with 
M2 = I.IMq, Pi = 1.04 d and 5 = 0.005 - 0.009. We also indicate in Fig. 3 whether 
the accretion disks in the LMXBs are thermally and viscously stable, according to the 
stability criterion for a mixed-composition (X = 0.1, Y = 0.9) disk from Lasota et al. 
(2008) We use the symbols x, *, and + on the evolutionary tracks to denote where 
the hydrogen composition X of the donor becomes 0.3, 0.2, and 0.1, respectively. 
The positions of UCXBs are marked in these two figures with circles and triangles 
for persistent and transient sources, respectively. Besides them, we also include 18 
NS LMXBs with known P and Mi (data are taken from Liu et al. 2007; Watts et al. 
2008; Heinke et al. 2009). 

A comparison between our CB disk-assisted binary models and the observations 
of (compact) NS LMXBs suggesting that it is possible to form UCXBs from normal 
LMXBs. We note that three of the UCXBs are in globular clusters, indicating low 
metallicities in these systems. However, from our calculations we find that change 
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Figure 2: Evolution of the orbital period for different values of the CB disk parameter 
5. 



of metallicities does not significantly affect the binary evolution when the CB disk is 
involved. 



5 Conclusion 

During mass transfer in LMXBs, a CB disk may be formed as a result of mass outflow 
from the accretion disk, and has been invoked as an efficient process for the removal 
of orbital angular momentum (Taam & Spruit 2001). We propose a scenario for the 
formation of UCXBs from L/IMXBs with the aid of a CB disk in this work. The 
suitable binary parameter space (M24 and Pi) with reasonable choice of the CB disk 
parameter S for the formation of UCXBs within 13.7 Gyr is found to be significantly 
larger than in previous "magnetic capture" model (van der Sluys et al. 2005). This 
difference is caused by the fact that the bifurcation period is considerably increased 
if the CB disk is included. 

This work was supported by Natural Science Foundation of China under grant 
10873008 and National Basic Research Program of China (973 Program 2009CB824800). 
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